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' We review cm and mm observations of thermal molecular line emission from high 

' redshift QSOs. These observations reveal the massive gas reservoirs (lO^" to 10" 

^ Mq) required to fuel star formation at high rates. We discuss evidence for ac- 

, five star formation in QSO host galaxies, and we show that these high redshift, 

FIR-luminous QSOs follow the non-linear trend of increasing Lj?//j/L'(CO) with 
increasing hpm. We conclude with a brief discussion of the recent CO detection 
of the most distant QSO at z = 6.42, and its implications for cosmic reionization. 



o 

^ ■ 1. Introduction 

Over the last few years, the study of high redshift QSOs has been revolu- 
. . tionized in three ways. First, wide field surveys have revealed lOO's of high z 

rS \ QSOs, right back to the epoch of cosmic reionization {z > 6; e.g., Fan et al. 

' 2003). Second, it has been shown that most (all?) low redshift spheroidal 

galaxies have central super-massive black holes (SMBH), and that the black 
hole mass correlates with bulge velocity dispersion. This Mb/j-cTv corre- 
lation suggests coeval formation of galaxies and SMBH, thereby making 
SMBHs a fundamental aspect of the galaxy formation process (Gebhardt 
et al. 2000). And third, mm surveys of high redshift QSOs find that 30% 
of the sources are 'hyper-luminous infrared galaxies' {LpiR — 10^'^ Lq), 
corresponding to thermal emission from warm dust, and that this fraction 
is independent of redshift out to z = 6.4 (Omont et al. this vol.). If the 
dust is heated by star formation, the implied star formation rates are ex- 
treme (> 10^ Mq year ), consistent with the formation of a large elliptical 
galaxy on a dynamical timescale of 10^ years. On the other hand, the FIR 
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luminosity constitutes typically only 10% of the bolometric luminosity of 
the sources, such that dust heating by the AGN remains an alternative. 
Demographic studies show that SMBHs acquire most of their mass during 
major accretion events marked by the QSO phenomenon (Yu & Tremain 
2002). 

Molecular line observations (typically CO) of FIR- luminous high z QSOs 
have revealed large gas masses in most cases observed to date (see Table 
1). Such gas reservoirs arc a prerequisite for star formation models for dust 
heating in FIR-luminous high z QSOs. The typical gas depletion timcscales 
are of order 10'' to 10^ years, if the dust is heated by star formation. In 
this review we consider this question in more detail. We restrict ourselves 
to z > 2 QSOs (sec also Barvainis 1999; and sec Scovillc ct al. 2004 
and Sanders & Mirabel 1996 for observations of lower redshift QSOs). We 
assume a standard concordance cosmology. 

2. Statistics 

Table 1 shows all the published detections of CO emission from sources at 
z>2 (only the most recent reference is given). Most of the sources are AGN 

(QSOs, radio galaxies = RG), since in these cases optical spectroscopic 
redshifts are available, although a number of optically brighter 'submm 
galaxies' (SMM) have now been detected in CO emission (Neri et al. 2003). 
About half the sources are strongly lensed (L). Most sources have been 
studied in the higher order transitions (> CO 3-2), although at z > 4 the 
lower order transitions become accessible to cm telescopes such as the VLA. 

Figure 1 shows the correlation between FIR luminosity and velocity 
integrated CO 1-0 luminosity {Pcoii-a) ^ ^™ '^^^ P*^^) '^'-"^ sources at low 
and high redshift (Beelen et al. 2004). The sources in Table 1 consitute 
most of the sources with Lpiji > 10^^ Lq in Figure 1. For sources without 
CO 1-0 measurements, the 1-0 luminosity was calculated assuming constant 
brightness temperature. For the high redshift sources the FIR luminosity 
is given approximately by: LpiR = 4 x lO^^(S'25o/mJy)L0, appropriate for 
a typical ULIRG SED (Omont et al. al. 2003; note that in Fig 1 the values 
were calculated using the measured multifrequcncy SEDs where available). 
For dust heating by star formation, the total star formation rate (SFR; from 
0.1 to 100 Mq) is given by; SFR = 4 x IQ-^^LpiR M© year-^ Gas masses 
can be derived from; M(H2) = X x -£'c7(9(i_o) j where X = 4 for typical spiral 
galaxies, and X = 0.8 for ULIRGs (Downcs & Solomon 1998). Note that 
X = 0.2 is the minimum (ie. optically thin) value, assuming solar C and O 
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J1909+722 


RG 


3 


534 


4- 


3 


1.6 


4 


Papadopoulos 2000 


4C60.07 


RG 


3 


788 


1- 





0.24 
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Carilli 2002 


B1202-0725 


QSO 
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Carilli 2002 
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abundances, and warm gas (70 K). 



The straight hnes correspond to hnear relations, with a ratio of 
Lfir/L'qq = 40 (dash) or 300 (sohd). The correlation is clearly non- 
linear, and is better described by a power-law of the form: LpiR oc LJ^q^-^. 
This non-linearity has been interpreted as an increasing 'star formation ef- 
ficiency' (= SFR/gas mass) with increasing SFR (Gao & Solomon 2003). 
An alternate hypothesis is that the AGN dominates dust heating at the 
highest luminosities. But this latter hypothesis begs the question: why 
would there be any correlation at all? Meaning that if two very different 
physical mechanisms are involved, one might expect a discontinuity in this 
relationship going from starbursts to AGN. 

3. Examples 

Figure 2 shows examples of the CO excitation for three z > 4 QSOs. The 
line strengths are consistent with constant brightness temperature (T^) up 
to CO 5-4. This is similar to what is seen for nuclear starburst regions, 
such as in M82, and implies dense (> 10"* cm~"^), warm (> 50 K) gas. The 
recent detection of HCN emission from the Cloverleaf quasar aX z = 2.56 
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Figure 1. The relationship between Lpiji and -^c'O(i-O) (^^^^^^ al. 2004). 

supports the idea of dense, warm (star forming?) molecular gas in these 
systems (Solomon et al. 2003). 



. B1202-0725 z=4.69 - circles ^Constant 
B1335-0417 z=4.41 - squares / 
Jl 148+5251 z=6.42 - triangles 




LVG(70K,2e4cm- 



I I 

B1202-0725 z=4.69 — circles 
B1335-0417 z-4.41 — squares 
J1148+5251 z=6.42 -- triangles 



200 400 600 800 1000 

Frequency (GHz) 




Figure 2. Left: The CO ladder for three high z QSOs. Right: The radio to IR SEDs. 
The solid line is the fit for M82. 
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Also shown in Figure 2 are the radio to IR SEDs of these sources. Again, 
the SEDs are similar to those seen in typical star forming galaxies. VLBA 
imaging of the 1.4GHz continuum emission from J1409+5628 at 2; = 2.58 
implies intrinsic brightness temperatures < 10^ K, consistent with the non- 
thermal emission expected from a star forming galaxy, but inconsistent with 
that expected from an AGN (Beelen et al. 2004). 




Figure 3. The optical (left) and CO 2-1 (right) images of the gravitationally lensed 
QSO J2322+1944 at z = 4.12 (CariUi et al. 2003). 

Figure 3 shows the optical and CO emission from the z = 4.12 QSO 
2322+1944. The source is strongly gravitationally lensed, appearing as a 
double QSO in the optical with 1" separation (Djorgovski et al. in prep), 
and a complete 'Einstein ring' in the CO 2-1 emission. A similar ring is seen 
for the 1.4 GHz continuum emission (Carilli et al. 2003). These results can 
be modeled in the source plane as a starburst disk surrounding the QSO 
with a radius of about 2 kpc. J2322+1944 represents perhaps the best 
example of a coeval starburst+AGN at high redshift. 

We consider briefly the question of why only 30% of high z QSOs are 
FIR luminous? Unfortunately, sensitivity limits of current (sub)mm surveys 
cannot rule out a continuum of FIR luminosities extending below a few 
mJy, although deep radio surveys suggest that this may not be the case 
(Petric et al. 2004 in prep). If there really are two types of QSOs (FIR- 
luminous and FIR-quiet), then the 30% fraction could represent a relative 
duty cycle for star formation vs. QSO activity. Alternatively, there may be 
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two types of major accretion events for SMBHs: those that have associated 
star formation, and those that do not. 

4. Cosmic Stromgren spheres 

The epoch of reionization (EoR) represents a fundamental benchmark in 
cosmic structure formation, corresponding to ionization of the neutral IGM 
by the first luminous sources. The recent discovery of Lya absorption by 
the neutral IGM toward z > 6 QSOs (the Gunn-Peterson effect; Fan et al. 
2003) implies that these sources are situated at the end of the EoR. 




40 y=Si a::: — i i 40 1- — h \ 1 l 

11 48 17.5 17.0 16.5 16.0 11 48 17.5 17.0 16.5 16.0 



RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Figure 4. CO 3-2 emission from tiie highest redshift QSO, J1148+5251 at z = 6.419 
(Walter et al. 2003). The left frame shows tlie velocity integrated line emission at 
1.5" resolution at 46.7 GHz. The right frame is an 'off' channel, showing the lack of 
continuum emission. The contour levels are: -0.1,0.1,0.2,0.3,0.4,0.5 mjy/beam. 

Bertoldi ct al. (2003) show that the 30% fraction of FIR-luminous 
QSOs remains constant into the EoR, including the most distant QSO 
known, J1148-I-5251, at 2: = 6.42, which has an Lpjji = 1.2 x 10^^ Lq. 
CO emission has also been detected from J1148-I-5251 using the VLA and 
the PdBI (Figure 4; Walter et al. 2003; Bertoldi et al. 2003), with an 
implied molecular gas mass of 2 x 10^" M©. The presence of a large mass 
of heavy elements and dust in a galaxy just 0.8 Gyrs from the big bang 
raises interesting issues for ISM enrichment, and suggests that active star 
formation started at 2; > 10 in the host galaxy of this QSO. Also required 
is a dust formation mechanism involving high mass stars and/or supernova 
remnants. 

A particularly interesting result for J1148+5251 is the difference be- 
tween the host galaxy redshift of 6.42 (as determined from the CO line) 
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and the on-set of the Gunn-Peterson absorption trough at z = 6.32. We 
note these CO observations are the only measurements relating to the 
host galaxy properties (eg. the exact redshift), as opposed to the AGN 
properties. This redshift difference implies that the QSO must be sur- 
rounded by a ionized sphere of physical radius R=4.7 Mpc, presumably 
formed by the radiation from the QSO itself. Hence, in J1148-I-5251 we 
are witnessing the process of cosmic reionization. This is a 'time bounded' 
Stromgren sphere, and implies a lifetime for the recent QSO activity of: 
tqso = 10^(R/4.5Mpc)^f(HI) years, where f(HI) is the (volume averaged) 
IGM neutral fraction (Walter et al. 2003; Wyithe & Loeb 2003a). 

Wyithe & Loeb (2003b) have used this result, plus models of QSO 
formation, to argue (statistically) that the neutral fraction of the IGM 
must be substantial at z = 6.4: f(HI) > 0.1, otherwise the QSO lifetimes 
become unreasonably short. This neutral fraction is much larger than the 
best current lower limits set by the GP effect of f(HI) > 0.001, and argues 
for 'fast reionization' at z ~ 6 (Gnedin 2000). The implication is that 
f(HI) changes from < lO"'' at z = 5.7 to > lO'^ at z = 6.4, ie. the 
neutral fraction of the universe changes by three orders of magnitude over 
a timescale of only 0.1 Gyr. 
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